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Abstract 

A new nine-dimensional (9D), ah initio electric dipole moment surface (DMS) of methane in 
its ground electronic state is presented. The DMS is computed using an explicitly correlated cou- 
pled cluster CCSD(T)-F12 method in conjunction with an F12-optimized correlation consistent 
basis set of the TZ-family. A symmetrized molecular bond representation is used to parame- 
terise the 9D DMS in terms of sixth-order polynomials. Vibrational transition moments as well 
as band intensities for a large number of IR-active vibrational bands of ^^CH4 are computed by 
vibrationally averaging the ah initio dipole moment components. The vibrational wavefunctions 
required for these averages are computed variationally using the program TROVE and a new 
'spectroscopic' ^^CH4 potential energy surface. The new DMS will be used to produce a hot 
line list for 12CH4. 
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1. Introduction 



Methane plays an important role in atmospheric and astrophysical chemistry. Its rotation- 
vibration spectrum is of key importance for models of the atmospheres of bodies ranging from 
Titan to cool stars. However the lack of precise data on methane spectra, particularly at higher 
teinperatures, has severely limited models for atmospheres as diverse as Jupiter l|, exoplanets 
[2,3] and brown dwarfs 4| . Any temperature-dependent model of the methane spectrum requires 
a reliable treatment of the associated transition intensities. The construction of a reliable electric 
dipole moment surface (DMS) which can be used to generate such intensities concerns us in this 
paper. 

Theoretically, a number of ab initio potential energy surfaces (PES) computed at different 
levels of theory have been reported [5 -10]. The surface by Schwenke and Partridge has 
been used in many benchmark studies as well as for accurate predictions of vibrational and 
ro-vibrational energies of CH4. However, to our knowledge no accurate and comprehensive ab 
initio DMSs of CH4 have been published so far. In some cases ab initio DMSs were used, but not 
actually provided, for example by Oyanagi et al. (at the CCSD(T)/cc-pVTZ and MP2/cc- 



pVTZ level of theory) and Warmbier et al. (RCCSD(T)/aug-cc-pVTZ). An empirical DMS 
of limited accuracy was developed by Hollenstein et al. [12l |. 

The goal of this paper is to bridge this gap and provide an accurate and detailed ab initio 
DMS of CH4. To this end we employed an explicitly correlated (F12) coupled cluster method 
in conjunction with an F12-optimized basis set to generate a nine-dimensional (9D) DMS of 
CH4. This surface is then used to compute vibrational transition moments and vibrational 
band intensities of ^^CH4 for a large number of IR-active transitions in the range between 
and 10000 cm~^. This is done by vibrationally averaging the DMS over vibrational functions 
of ^^CH4 computed through variationally solving the nuclear motion Schrodinger equation with 
the TROVE program {l^. This program has been used successfully for accurate and extensive 
ro-vibrational calculations on different tetratomic molecules 14|-|2ll|. The present work is the 
first application of TROVE to a larger system. Other variational approaches capable of treatin; 
this pentatomic molecule include the Lanczos-method developed by Carrington and Wang 



251], MULTIMODE ^^\2^, and GENIUSH [sdl as weh as the methods of Halonen 



and Tennyson 
employed by Yu 



3211 . Schwenke 331 ] . and Nikitin 



35, 



22 



3l|], Xie 



341 ]. The rigorous quantum dynamics algorithm 



361 ] should be also mentioned as an alternative accurate approach employed 



for calculating vibrational energies of CH4. 

On the experimental side, electric transition dipole moments for CH4 are the subject of a 
large number of publications 



37 



471 ] ■ where effective dipole moment parameters were derived 
2 



either through fits to intensity measurements or using perturbation theory. The most recent 
and comprehensive collection of the effective dipole moment parameters describing excitations 
in the region up to 4800 cm~^ was reported recently by Albert et al. Many of these studies 
also provide band intensities. We use some of these data to validate our vibrational transition 
moments as well as our DMS, where possible. 

The paper is structured as follows. In Section [2] the ab initio method used for generating 
DMS of CH4 is described and an analytical representation for this DMS is introduced. The 
PES is detailed in Section [3l The variational approach is described in Section 2] and theoretical 
vibrational transition moments and band intensities are reported in Section [5j Finally, Section [6] 
offers some discussion and future prospects. 



2. Ab initio calculations 

We used the recently proposed explicitly correlated F12 coupled cluster method with sin- 
gle and double excitations including a perturbational estimate of connected triple excitations, 

1 the corresponding F12-optimized correlation-consistent 



CCSD(T)-F12c [48!l in conjunction wit: 
basis sets, namely aug-cc-pVTZ-F12 49(] (augmented correlation-consistent polarized valence 
triple-zeta basis) as implemented in Molpro 2010 jsol . 

Methane has nine vibrational degrees of freedom and thus requires a very extensive grid 
of geometries in the ab initio calculations. We have computed three ab initio DMSs (one for 
each Cartesian component) on a fine grid of 114,000 geometries covering the energy range up to 
50,000 cm~^ with C~H distances ranging from 0.75 to 2.0 A, and interbond angles ranging from 
50 to 140°. The finite field method, with a field of 0.005 a.u., was used to compute the dipole 
moment component as first derivatives with respect to the external electric field in the dipole 
approximation. 

In order to represent the ab initio dipole moment vector fl of CH4 analytically we use 



the so-called symmetrized molecular bond (SMB) representation 



where /z is projected onto 



three symmetrized combinations of vectors associated with the four molecular bonds C-Hj (i 



1, 2, 3, 4). An analogous approach was used to represent the ab initio dipole moment of NH3 [51| 
(see also Ref. [s] for a non-symmetrized MB representation of the methane DMS). The three 
components of the dipole moment of CH4 span three components of the F2 representation, i.e. 
F2x, F2y, and in the irreducible representation of the Td(M) molecular symmetry group |52l ]. 



We define the following three symmetrically independent reference vectors iir spanning the same 
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representation F2: 

2 \\ri\ \r2\ IrsI \r4\J 

2 \|ri| \r2\ Iral \r4\J 

^ ( n f2 n f4 \ 



where is a molecular bond vector pointing from C to Hj [i = 1,2,3,4). The dipole moment 
vector jl can now be represented as 

/I = n^a, ^^F2^ + np^y fJ-Fiy + /^Faz , (4) 

where the projections f^F2a i'^ — ^lUi^) ^-^e the dipole moment functions also spanning the F2 
representation. Each component fJ.F2a is then expanded as 



ll,J2v:«9 

in terms of the following set of nine independent internal coordinates: 



= (r,-re)e-^('^-^<=)', i = l,2,3,4, (6) 

^5 = ^=(2ai2 - ai3 - ai4 - "23 - "24 + 2a34), (7) 
V 12 



?6 = ^(ai3 - ai4 - "23 + "24), (8) 

?7 = ■^(«24-ai3), (9) 

6 = ;^("23 - "m), (10) 

?9 = "^("34 - 012), (11) 



where rj are the bond lengths and Uij are the interbond angles. Following Hollenstein et al. [JJ 
we have introduced the factor exp (^—^{ri — Tg)^) in order to keep the expansion in Eq. ([5]) from 
diverging at large r^. Here is a reference expansion center which is conveniently taken at a 
value close to the molecular equilibrium. 

The properties of the components /iFa^ associated with F2 impose symmetry relations be- 



tween the expansion coefficients /if^ j^. 
standard symmetry transformation rules 



hese relations can be reconstructed by applying the 



52 



53I ] between the F2x, F2y, F2z components of F2. 



We employed the Maple program function 'solve' to find a set of independent coefficients ^f^ 



■,«9 



54( 1 for details on the com- 



by solving an over-determined system of linear equations (see Ref. 
putational procedure). For example, to first order the dipole moment functions are given by 

4 



the expansions 

= M7 + m(6-e2 + 6-C4) + o(ei)..., (12) 

f^F.y = M8 + w(ei-e2-6 + C4) + o(ei)..., (13) 

t^F2. = M9 + m(ei + e2-6-^4) + o(ei)..., (m) 

where the shorthand notations ^7 = Aio,o,o,o,o,o,i,o,o /^i ~ /"i,o,o,o,o,o,o,o,o used. The 
complete expansion set is given in the supplementary material to this paper. A full sixth-order 
symmetrized expansion requires 680 symmetrically independent dipole moment parameters in 
total. Taking into account the symmetry relations, the expansion Eq. ([5]) can be written as: 

21,12, .■■1*9 

where {^^^^2^ . . -Cg^}^^" is a symmetrized combination of different permutations of the internal 
coordinates consistent with T(j(M) group operations and spanning the symmetry. 

The actual values of expansion parameters /x"^ 13 ig corresponding to our ab initio DMS 
were determined through least-squares fittings to the 114,000 ab initio dipole moment values. 
Some parameters with large standard errors that exhibit strong correlation were constrained to 
zero. The ab initio data in these fits were weighted using the following expression proposed by 
Partridge and Schwenke 55l |: 



tanh -0.0004 cm x (Vi - 18 000 cm^^) + 1.00000002 , , 

Wi = N ^ ^— '-^ , 16 

2.00000002 ' ^ ^ 

where is a normalization factor defined by ^ ■ wi = 1 and Vi is the corresponding ab initio 
energy at the ith geometry (in cm~^), measured relative to the equilibrium energy and com- 
puted using the same level of theory. The ab initio energy Vi is weighted by the factor Wi in 
the PES fitting; these weight factors favour the energies below 18 000 cm~^. In the end 296 
parameters were needed to obtain a weighted (wi) root-mean-squares (rms) error of 0.00016 D. 
The expansion parameters are given in the supplementary materials of the paper along with a 
fortran 95 routine incorporating the corresponding analytical form. 

3. Potential energy surface 

In this work we use a new 'spectroscopic' PES generated as follows. First, an ab initio 
PES was gen erated on the same grid of 114,000 geometries employing the same CCSD(T)-F12c 



approach 



48( 1 as for the DMS, but with the larger aug-cc-pVQZ-F12 basis set 



491 ]. The frozen 



core approximation was applied in the coupled cluster calculations. We did not include core- 
valence corrections and higher-order coupled cluster corrections because they cancel to a large 

5 



extent (see the work by Yachmenev et al. 



18|), and also because this would be very costly 



with the resources available to us. Relativistic corrections were computed using the Douglas- 
Kroll method as implemented in Molpro 2010. With this method we obtained the equilibrium 
bond length = 1.08734 A, which is close to the value = 1.08595 A derived in a combined 



experimental/ a& initio analysis 



5a]. The ah initio potential energies have been represented using 



a symmetrized analytical expansion in terms of the four stretching coordinates 1 — exp(— a(rj — 
Te)) {i = 1,2,3,4) and five bending coordinates (j = 5,6,7,8,9) from Eqs. ([6 HTT]) . The 
symmetry relations between potential parameters in the Td(M) group were derived employing 
Maple with the procedure described above. In total 287 symmetrically independent parameters 
are needed for a full sixth-order expansion. These expansion parameters were obtained from a 
fitting procedure using the weighting scheme given in Eq. ()16p . Some parameters showing large 
standard errors were constrained to zero, which resulted in 268 potential parameters needed to 
get a weighted rms error of 0.3 cm~^. 

To improve the quality of the ah initio PES, we slightly refined it by fitting to the 'experi- 



57l | database for J = 0, 1, 2, 3, 4. 



mental' energies of ^2CH4 extracted from the HITRAN 2008 

In these fittings the variational program TROVE [3] was used, closely following the approach 
described in detail by Yurchenko et al. 5^. This new 'spectroscopic' PES for ^^CH4 will be 
the subject of further improvements and a future publication. Both the initial ah initio and re- 
fined potential parameters are given as supplementary material along with a fortran 95 program 
containing the corresponding potential energy function. 



4. Variational calculations 

The vibrational wavefunctions and energies of ^^CIl4 were calculated using the variational 
program TROVE In this approach the Hamiltonian operator is represented as an expansion 
around a reference geometry taken presently at the molecular equilibrium. The kinetic energy 
operator is expanded in terms of the nine linearized coordinates ^f, which are linearized versions 
of the internal coordinates {i = 1,...,9) defined above; for the potential energy function 
we use 1 — exp(— a^f) {i = 1,2,3,4) for the four stretching and .^j [j = 5,6,7,8,9) for the 
five bending linearized coordinates. The latter expansions were applied to the refined PES 
introduced above. The kinetic and potential energy parts were truncated after the sixth and 
eighth order, respectively. 

To construct the basis set, TROVE employs a multi-step contraction scheme based on the 
following polyad truncation. The polyad number in case of CII4 is given by 

P = 2{vi +V2 + VZ+ Vi) +V^ + V^ + Vj + V^+ Vg, (17) 
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where Vi is a vibrational quantum number associated with a one-dimensional primitive basis 
function (j^vii^f) {i = 1, 2, 3, . . . , 9). At Step 1, nine sets of (/>j;,('?f) generated using the 



Numerov-Cooley 



59 



60] method by solving the ID Schrodinger equations for each of the nine 
modes separately. The corresponding ID Hamiltonian operators are obtained from the 9D 
Hamiltonian operator (J = 0) by freezing all but one vibrational coordinates at their equilib- 
rium values. At Step 2 the 9D coordinate space is divided into the three reduced sub-spaces, 

(i) ^2' ^ll' (ii) ?5i'?6' ^^^'^ (iii) ■^yjCljCg- This division is dictated by symmetry: each of 
the sub-spaces is symmetrically independent and thus can be processed separately. For each of 
these sub-spaces, the Schrodinger equations are solved for the corresponding reduced Hamilto- 
nian operators employing as basis the products of the corresponding primitive functions 4>Vi{^i) 
generated at Step 1. The sub-space bases are truncated using the condition P < Pmax, where 
P is given by Eq. p!7|) . In this work we use Pmax = 12, i.e. our reduced sub-space basis sets are 
defined by 

"f^l + ^^2 + ^3 + < 6, (18) 

V5 + ve < 12, (19) 

V7 + V8 + Vg < 12. (20) 

The three sets of eigenfunctions ^-JJ (Cf, d> ^1, d), ^S^Cei^D, and (^^,^1,?^) resulting from 
these solutions are symmetrized and assigned local mode quantum numbers (i) vi,V2,V3,V4, 

(ii) V5,vq, and (iii) ^7,^8,^9 and vibrational symmetry T^n,. The final vibrational basis set is 
formed from products ^^'^^ x ^^"^ x which are contracted with P < 12 and symmetrized 
again using the standard technique for transformations to irreducible representations (see, for 



example, Ref. [5^). Thus our contracted vibrational basis set consists of ten symmetrically 
independent groups, one for each irreducible representation of the Td{M) group Ai, A2, Ea, 
Eb, Fix, Fly, Fiz, F2x, F2y, F2z- From the degenerate symmetries E,Fi, and F2, only the first 
component is processed. Therefore only five vibrational Hamiltonian matrices (one for each 
symmetry) are constructed and diagonalized. In the diagonalizations we use the eigensolver 
DSYEV from LAPACK as implemented in the MKL libraries. The dimensions of the matrices 
to be diagonalized are 2190, 1725, 3901, 5481, and 5940 for Ai, A2, E, Fi, and F2, respectively. 
The results are collected in Tables [TH3l where the calculated energies are also compared with the 



available experimental values from 



61 



-63||. It should be noted that our basis set with P < 12 

does not provide fully converged energies. The same comment applies to our new 'spectroscopic' 

PES: because of the limitations of the method, this PES can only guarantee the energies given 

in Table [1] when used with the same approach and basis set employed in the refinements, i.e. 

7 



-Pmax = 12- However the vibrational transition moments reported in the next section appear to 
be rather insensitive to the size of the basis set and the shape of the potential. 



5. Transition moments 

The transition moment connecting two states i ('initial') and / ('final') can be defined as 
the vibrational average of a dipole moment component fia {ct = x, y, z): 



r' r" 



(21) 



where fia (" = x, y, z) are the Cartesian components of the 



dipole moment in the molecule- 

r' r" 
64j): ^ • " and ^'r'' are vibrational 



fixed frame (defined in TROVE by the Eckart conditions 
eigenfunctions obtained variationally, and Tg represents a component of one of the ten irreducible 
representations Ai, A2, Ea, E^, Fix, Fiy, Fi^, F2x, F2y, -F22 of rd(M). It should be noted 
that the transition moments defined by Eq. (j2ip depend on the choice of the molecular-fixed 
coordinate system [65(] as well as on the particular choice of the transformational properties of 



the irreducible representations E,Fi, and F2 (we use irrep-matrices from Hougen's monograph 
66t|). The following quantity, averaged over the Cartesian components a = x,y,z and any 
degenerate components of the wavefunctions, is independent of these choices: 



fl^f = /i(^'r';/,r") 



(i,r';/,r") 



+ 



-(i,r';f,r") 



+ 



■-,(i,r';/,r") 



(22) 



where 



1^. 



(i,r';/,r") 



(23) 



a,b 



(i,r';/,r") 

X ) 



and the vibrational transition moment /x*-^ is defined as the length of the vector [/i, 
(i,r';/,r") (i,r';/,r")n 

A number of calculated values of the vibrational transition moments /i*-^ as well as matrix 
elements Jii''^ ''^'^ ^ of ^^CH4 corresponding to transitions from the three lowest vibrational states 



('^Oi ^2, and Vi) are listed in Tables [TH31 The non-zero matrix elements /X; 
for each transition satisfy the following set of relations: 



(i,r';/T") 



in Tables [IH31 



H-y 

2/2f-^^") : 
l-l'X 



_r-,iAv,F2y) 
h^Z ) 

,-.(i^2x;F2.) 
H'Z , 

^ (Fa, 
-{F2y;Fi^) 



-iF2^;Fi,)_ 

H'X 

-(Fi^;Eb) _ 2 r-S^iy;Ea)_ r,-(Fiy;Et) __ 2 -iFi.;Ea)_ 

rX ~ y/3^y ~ ~ 



2/2f'^'^"^ 

r,iF2y;Fi,) 
H-y 



_J_-{F2,;E,) 

. -(i^2.;Fi,) 
- 

= 2/21^-'^"), 



(24) 



{F2z;Fiy) 



where indices i,f are omitted for simplicity. Therefore Tables [TH3] list only one value for each 
pair of states. A complete list of 47861 vibrational transition moments and matrix elements 



Jla ' of ^^CH4 is provided as supplementary material for all allowed transitions in the 
frequency range from to 10000 cm^^ with lower state energies below 5000 cm"^. 

Because of the high symmetry of methane many vibrational transition moments p}^ vanish. 
For example, fi^^ = between Ai eigenfunctions of the ground vibrational state (GS) and any 
vibrational states of Ai, A2, E, and Fi symmetry. This includes transitions within the GS and 
the transitions giving rise to the fundamental bands uii^Ai) and U2{E). All these transitions 
have zero vibrational transition moments and are thus dipole forbidden. Certainly this does not 
rule out ro-vibrational transitions within these bands although such transitions are often also 



referred to as 'forbidden'. At a formal leve^ 
with the correct symmetries are allowed 67 



any transitions connecting rotation- vibration states 



68( 1 . In practice, centrifugal distortions arising from 



rotation-vibration interactions can introduce i<2-vibrational contributions into, for example, the 

GS rotational wavefunctions, thus producing non-zero ro-vibrational line strengths. Conversely, 

'allowed' ro-vibrational transitions are based only on non-zero vibrational transition moments 
(i r' ■/ r") 

/2a' ' , such as the ones given in Tables [TH3l 

The transition moments presented here can be used to derive the effective dipole moment pa- 



rameters of CH4, which are important for modelling intensities of this molecule (see Ref. 



47| 



for the most recent collection of the effective dipole moment parameters). Examples of such 



an exercise can be found in Refs. [4 



69 1 . For example, the parameter 



0(0,0A1) 

OOOOA1-0001F2 



0.098677(88) D [47l is in good agreement with our value VS/ii^ " '^''"^'^'^ = 0.098314 D. Similarly, 



0(0,0A1) 



OOOOA1-OOIOF2 



0.094622(47) D [47|] agrees weU with our value 



■,.s.,Av,vz,F2) 



0.093699 D 



(see Tabled]). Here \/3 is a conversion factor associated with the Dyad-GS and Pentad-GS zero- 



38|, 



order effective dipole constants in the tensorial representation (see, for example, Refs. 
For higher vibrational order parameters the correspondence between our transition moments 
and effective dipole moment parameters is less straightforward owing to the complexity of the 
effective dipole moment operators which are usually given in a tensorial representation. 

From the transition moments computed we also generated vibrational band intensities as 
given by Brown et al. [38l | and 39[: 

o_3in— 36 T 



3/icQ^ T " y -*^' ^^^^ 

where and E\ are the vibrational transition frequency and the lower state energy, respectively, 
Lochschmidt's number L = 2.686754 x 10^^ cm-^, Tq = 273.15 K, k is the Boltzmann constant, 
and Qt, = 1.002 is the vibrational partition function for T = 278.15 K. Some of the computed 



values of S'vib are also illustrated in 
band intensities from Refs 



39|, 



3S 



42 



'ab: 



es [TH51 where they are compared to the experimental 



45144 7l|. where available. The agreement is generally very 



good. In some cases, particularly for combinational bands, there is some ambiguity over how 
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to distribute the individual lines between transitions with the same normal-mode vibrational 
quantum numbers (see also footnotes to Tables [TH3|) • This introduces an extra uncertainty into 
these comparisons. 

As an illustration of a possible application of our transition moments, Fig. [1] shows the 
vibrational band intensities calculated at T = 298.15 K (empty blue squares) and T = 1500 K 
(filled red circles) using Eq. ([25]) {Qy = 3.038 at T = 1500 K). 

6. Conclusion 

In this work we present a new ab initio DMS for CH4 which is used to generate the vibrational 
transition moments and band intensities of ^^CH4 for a large number of vibrationally allowed 
transitions. These values should be helpful for predicting intensities for individual bands. The 
dipole matrix elements can also be used for predicting or evaluating effective dipole moment 
parameters. The analytical representations for the electric dipole moment and potential energy 
functions of CH4 developed in this work have correct symmetry properties of the T^{M) group 
and should be useful for representing similar objects obtained using diff'erent methods or levels 
of theory. 

Comparisons of calculated and experimental (or experimentally derived) band intensities 
demonstrate the good quality of our new DMS of CH4. We are planning to use this DMS in 
calculations of a high temperature line list for methane within the framework of the ExoMol 
project [7^ (see www.exomol.com), which has identified methane as a key target species. Such a 
line list will be important for modelling molecular opacity in atmospheres of (exo-)planets and 
cool stars. 
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Table 1: Vibrational transition moments, fl'-^ in D, individual matrix elements /i^'"'^^'''^'^^^-' (D), and band in- 
tensities, Sif in cm^^atm^^, for ^^CH4 for transitions from the vibrationally ground state. The calculated and 
experimentally derived term values Ef of the upper states are given in cm^^. 



r 


State-f 


pobs 
f 


pcalc 
f 




-if 


cr(calc} 
if 


if 


F2 


U 


U 


U 


1 


1310.76 


1310.87 


-0.05676 


0.09831 


129.323 


127.68 


F2 











2 


2614.26 


2614.32 


0.00408 


0.00706 


1.334 


1.05 




U 


1 
i 


U 


1 


2830.32 


2830.28 


0.01038 


0.01797 


9.346 


6.63" 


F2 


U 


U 


1 


U 


3019.49 


3019.49 


-0.05410 


0.09370 


271.062 


269.92 


F2 











3 


3870.49 


3870.49 


-0.00231 


0.00401 


0.636 


59-1 


Fi 











3 


3930.92 


3930.81 


0.00166 


0.00288 


0.332 


O.14J 


Fo 


U 


1 


U 


r> 
Z 


4142.86 


4142.99 


-0.00171 


0.00295 


0.370 


0.34 


F2 


1 








1 


4223.46 


4223.52 


0.00775 


0.01343 


7.791 


7.84 


F2 








1 


1 


4319.21 


4319.43 


-0.00919 


0.01593 


11.201 


10.24 


Fo 
^ z 


U 





U 


1 


4348.72 


4348.97 


-0.00145 


0.00252 


0.281 


0.55-1 


F2 





2 





1 


4378.95 


4379.15 


-0.00013 


0.00023 


0.002 


O.23J 


F2 





1 


1 





4543.76 


4543.90 


0.00392 


0.00678 


2.138 


1.26'' 


Fo 
^ z 


U 


U 


U 


4 


5143.24 


5143.25 


0.00035 


0.00060 


0.019 


0.01'= 


F2 











4 


5211.29 


5210.65 


-0.00021 


0.00037 


0.007 




F2 


n 
u 


^ 


n 
u 


Q 


5370.52 


5370.37 


-0.00014 


0.00023 


0.003 




F2 





1 





3 


5429.58 


5429.07 


-0.00058 


0.00101 


0.056 




F2 





1 





3 


5445.12 


5444.87 


0.00011 


0.00020 


0.002 




F2 


1 








2 




5519.47 


0.00011 


0.00019 


0.002 




F2 








1 


2 


5588.03 


5586.76 


0.00156 


0.00270 


0.416 




F2 








1 


2 




5620.63 


-0.00074 


0.00129 


0.095 




F2 








1 


2 




5633.69 


-0.00102 


0.00176 


0.179 




F2 





2 





2 


5643.45 


5643.64 


0.00005 


0.00009 


0.000 




F2 





2 





2 


5668.60 


5668.70 


-0.00027 


0.00047 


0.013 




F2 


1 


1 





1 




5727.16 


0.00016 


0.00028 


0.005 




F2 





1 


1 


1 


5823.10 


5823.03 


-0.00213 


0.00368 


0.808 




F2 





1 


1 


1 


5844.00 


5843.72 


-0.00052 


0.00089 


0.048 




F2 


1 





1 







5861.92 


-0.00031 


0.00054 


0.018 




F2 





3 





1 


5867.66 


5868.23 


-0.00026 


0.00046 


0.012 




F2 





3 





1 


5894.12 


5894.54 


-0.00007 


0.00013 


0.001 




F2 








2 





6004.69 


6004.45 


-0.00301 


0.00522 


1.671 


1.63= 


F2 





2 


1 





6054.64 


6054.67 


-0.00055 


0.00094 


0.055 




F2 





2 


1 





6065.32 


6065.26 


-0.00055 


0.00095 


0.056 





" The total intensity of the 1/2 + 1^4 band, 9.94 cm^^atm^^ Albert et al. |47|| . which mcludcs both Fi and F2 sub-components 
agrees well with our Sgjg'jj^^ value. 

The total intensity of the 3i^4 band, 0.83 cm~latm~^ Albert et al. [47I . which includes Ai and Fi sub-components in 
addition to two F2 ones can be compared to the sum of the two 'calc' band intensities Sgg^^p^ . 

The total intensity of the 2u2 + Vi band which includes three sub-bands (Fi, F2 and F2) by Albert et al. [47I is 0.93 
cm^^atm^^. Hilico et al. [45I reported somewhat smaller value of 0.63 cm^^atm^^, which suggests a better agreement 
with our value for the F2 sub-band only. 

'' Compare to the estimate by Hilico et al. [ssl . 1.84 cm^^atm"'^, which is closer to our value. 
*= As estimated by Robert et al. 



11 



Table 2: Vibrational transition moments, in D, individual matrix elements /ij' °' ' ^'^ (a = x,y, z oi a,b 
and P — x,y,z) in D, and band intensities, 5*;/ in cm^^atm^^, for ^^CH4 for the transitions from the ^4 state. 
The calculated and experimentally derived term values Ef of the upper states are given in cm^^. 





Statc-f 




^calc 


_(i,rQ;/,F2f)) 


a 


13 




7 f 

' calc 


ry (calc) 


S)° ' [42] 


E 





1 








1533.33 


1533.39 


-0.00406 


a 


z 


0, 


,00704 


0.00013 




Ai 











2 


2587.04 


2586.98 


0.05046 




z 


0, 


,08741 


0.17802 


1 


F2 











2 


2614.26 


2614.32 


0.05718 


y 


X 


0, 


,14006 


0.46701 


>0.931'^ 


E 











2 


2624.62 


2624.61 


0.06530 


a 


z 


0, 


,11310 


0.30697 


J 


F2 





1 





1 


2830.32 


2830.28 


0.00641 


y 


X 


0, 


,01569 


0.00684 


0.00872 


Fi 





1 





1 


2846.07 


2846.06 


0.00007 


z 


X 


0, 


,00016 


0.000001 




Ai 


1 











2916.48 


2916.54 


0.00236 




z 


0, 


,00408 


0.00049 


0.000689 


F2 








1 





3019.49 


3019.49 


-0.01675 


a 


z 


0, 


,04104 


0.05265 


0.0481 


Ai 





2 








3063.65 


3063.61 


-0.00088 




z 


0, 


,00153 


0.00007 


0.0005 


E 





2 








3065.14 


3065.12 


0.00144 


a 


z 


0, 


,00250 


0.00020 




F2 











3 


3870.49 


3870.49 


-0.00383 


a 


z 


0, 


,00938 


0.00412 




Ai 











3 


Qnnn on 
o9U9.Z0 


3909.68 


n nn /i o /i 
-0.00424 




z 


0. 


nnvQ A 
,00/34 


O.OOzob 




Fi 











3 


3920.51 


3920.15 


-0.00320 


z 


X 


0, 


,00784 


0.00293 




F2 











3 


3930.92 


3930.81 


-0.00299 


a 


z 


0, 


,00732 


0.00257 




E 





1 





2 


4101.39 


4101.49 


-0.00973 


a 


z 


0, 


,01685 


0.01451 




Fi 





1 





2 


4128.76 


4128.73 


0.00890 


z 


X 


0, 


,02180 


0.02451 




Ai 





1 





2 


4132.86 


4132.83 


0.00923 




z 


0, 


,01599 


0.01320 




F2 





1 





2 


4142.86 


4142.99 


-0.00537 


a 


z 


0, 


,01315 


0.00897 




E 





1 





2 


4151.21 


4151.28 


0.00525 


a 


z 


0, 


,00910 


0.00430 




F2 


1 








1 


4223.46 


4223.52 


0.00244 


a 


z 


0, 


,00597 


0.00190 




F2 








1 


1 


4319.21 


4319.43 


-0.03798 


a 


z 


0, 


,09303 


0.47649 




E 








1 


1 


4322.18 


4322.21 


-0.04272 


a 


z 


0, 


,07399 


0.30170 




Ai 








1 


1 


4322.70 


4322.40 


-0.03326 




z 


0, 


,05761 


0.18292 




Fi 








1 


1 


4322.59 


4322.66 


-0.03760 


z 


X 


0, 


,09211 


0.46758 




F2 





2 





1 


4348.72 


4348.97 


-0.00451 


a 


z 


0, 


,01104 


0.00678 




Fi 





2 





1 


4363.61 


4363.61 


0.00076 


z 


X 


0, 


,00186 


0.00019 




F2 





2 





1 


4378.95 


4379.15 


-0.00057 


a 


z 


0, 


,00140 


0.00011 




E 


1 


1 








4435.12 


4435.18 


-0.00095 


a 


z 


0, 


,00165 


0.00015 




Fi 





1 


1 





4537.55 


4537.75 


0.00143 


z 


X 


0, 


,00351 


0.00073 




F2 





1 


1 





4543.76 


4543.90 


-0.00121 


a 


z 


0, 


,00297 


0.00052 




E 





3 








4592.03 


4592.12 


-0.00030 


a 


z 


0, 


,00052 


0.00002 




Ai 





3 








4595.51 


4595.66 


0.00000 




z 


0, 


,00000 


0.00000 





" This 'experimental' vibrational intensity of the 2u4 band [4^ correlates with our sj^'^2Ai ~^ '-'0002F2 '-'0002B ~ 0.952 
cm~^atm~'^. 
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Table 3: Vibrational transition moments, fl^^ in D, individual matrix elements p.^' °" ' {a = x,y,z and 
(3 = a,b) in D. and band intensities, Sif in cm^^atm"^, for ^^CH4 for the transitions from the V2 state. The 
calculated and experimentally derived term values Ei of the upper states are given in cm~^. 



r f 


State-f 


TT'obs 


IT' c ale 






p 


-if 
^calc 




o(calc) 


C.(obs) r 

i>if [42] 


F2 











2 


ZD14.ZD 
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r\ r\r\ 1 1 r\ 

0.00110 


z 


a 


U.0U191 


0.00002 


0.00008 


77 

-t<2 





1 





1 


Z53U.oi 


on 00 
2o3U.2o 


-0.UO/4U 


z 


a 


0.09942 


0.08000 




























I0.I45" 
J 


Fi 





1 





1 


2846.07 


2846.06 


-0.05718 


X 


b 
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0.08037 


F2 








1 





3019.49 


3019.49 


-0.01835 


z 


a 


0.03178 
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F2 











3 
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z 


a 
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Fi 











3 
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X 
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F2 
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z 


a 
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X 


b 
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0, 
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z 


a 


0.00507 


0, 
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F2 


1 








1 
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z 


a 
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0, 
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1 


1 
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z 


a 
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0, 
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1 


1 
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X 


b 
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0, 


,00025 
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z 


a 
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0, 
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X 


b 
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0, 
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1 
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z 


a 
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0, 
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4 
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z 


a 
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0, 
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F2 











4 
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z 


a 
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0, 
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4 
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X 


b 
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0, 
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3 
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-0.00262 


z 


a 


0.00454 


0, 
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b 
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F2 





1 





3 


5445.12 


5444.87 


-0.00064 
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" This 'experimental' vibrational intensity of the l>2 + 2^4 band |42l correlates with our 
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Figure 1: Vibrational band intensities calculated at T = 298.15 K (filled blue circles) and T = 1500 K (empty 
red squares) using Eq. 
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